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CO spectral line formation in the Sun has long been a source of consternation for solar physicists, as have the elemental abundances it seems to 
imply. We modelled solar CO line formation using a realistic, ab initio, time-dependent 3D radiative-hydrodynamic model atmosphere. Results 
were compared with space-based observations from the ATMOS space shuttle experiment. We employed weak 12 C 16 0, 13 C 16 and 12 C ls O lines 
\ from the fundamental (Av = 1) and first overtone (Av = 2) bands to determine the solar carbon abundance, as well as the 12 C/ 13 C and 16 0/ 18 

■ isotopic ratios. A weighted solar carbon abundance of log 6c = 8.39 ± 0.05 was found. We note with satisfaction that the derived abundance 

i-^ is identical to our recent 3D determination based on Ci, [Ci], C-> and CH lines, increasing our confidence in the accuracy of both results. 
Identical calculations were carried out using ID models, but only the 3D model was able to produce abundance agreement between different 
Q ■ CO lines and the other atomic and molecular diagnostics. Solar 12 C/ 13 C and 16 0/ ls O ratios were measured as 86.8+ 3 ' (6 >3 C = 30^) and 479^ 
' (6 ls O = 41**g), respectively. These values may require current theories of solar system formation, such as the CO self-shielding hypothesis, to 
(/} , be revised. Excellent agreement was seen between observed and predicted weak CO line shapes, without invoking micro- or macroturbulence. 

Agreement breaks down for the strongest CO lines however, which are formed in very high atmospheric layers. Whilst the line asymmetries 
^ (bisectors) were reasonably well reproduced, line strengths predicted on the basis of C and O abundances from other diagnostics were weaker 
• <-H , than observed. The simplest explanation is that temperatures are overestimated in the highest layers of the 3D simulation. Thus, our analysis 
' supports the presence of a COmosphere above the traditional photospheric temperature minimum, with an average temperature of less than 
^ , 4000K. The shortcoming of the model atmosphere is not surprising, given that it was never intended to properly describe such high layers. 
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1. Introduction 

1.1. CO in the Solar Spectrum 

Thousands of vibration-rotation lines from the fundamental 
1 2 + electronic band of CO are present in the infrared solar spec- 
trum obtained from space in the Atmospheric Trace Molecule 
Spectroscopy (ATMOS) space shuttle mission (see Sects. [3]and 
E). From 1350 to 2328 cnr 1 (i.e. 4.3 to 7.4 yum; fundamental 
bands from 1-0 to 20-19) and 3410 to 4360 cm- 1 (i.e. 2.3 
to 2.9 /mi; first overtone bands from 2-0 to 14-12), the so- 
lar spectrum looks more like a pure CO absorption spectrum, 
displaying not only 12 C ls O lines but also I3 C ie O, 12 C I8 and 
12 C n O isotopomeric lines. Thanks to their very high excita- 
tions (higher than in any laboratory spectra where rotational 
excitations are concerned), the solar CO lines have been used 
to derive a new set of highly accurate molecular constants 
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jFarrena et alJll99ltlSauval et alJll992l) . As the most sensitive 
indicator of solar photospheric temperatures, and formed over 
an extremely wide domain of optical depths (Fig. analy- 
sis of the CO lines offers a unique opportunity to test physi- 
cal conditions up to very high atmospheric layers, derive iso- 
topic abundan ce ratios and refine the solar carbon abundanc e 
JGrevesse&Sauva]|ll99ll 1 19921 1 1994 iGrevesse et allll995h . 
Only a limited number of clean CO lines are available in 
ground-based spectra however, with the Av = 1 sequence heav- 
ily masked by strong, broad telluric absorption and the Av = 2 
sequence also strongly polluted by terrestrial lines. 

1.2. CO Analyses and Problems 

Since the seminal observations of iNoves & HafJ il972h re- 
vealed limb CO core brightness temperatures below the clas- 
sical photospheric temperature minimum of 4500K, debate has 
raged over what CO absorption features in the solar spectrum 
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really tell us. lAvres & Testermanl J 198 lh found that the most 
natural explanation was an absence of any photospheric tem- 
perature minimum, requiring an outwardly decreasing temper- 
ature structure from the base of the photosphere right through 
the chromosphere. However, given that Can and Mgn line 
cores indicate c hromospheric temperatures we l l abov e the clas- 
sical minimum (lAvres & Lin skv 1976), Avres ( 1981) proposed 
the existence of a cool CO structure in the low chromosphere, 
coexistent with other areas of hotter gas responsible for the Ca n 
and Mg n core emission: a 'thermally bifurcated' low chromo- 
sphere; 

lAvres & Wiedemann! (|l989) went on to test the LTE as- 
sumption of the bifurcation model, finding that non-LTE ef- 
fects had <2% effect upon CO line cores, even as near to the 
solar limb as u — 0.1. Whilst these findings were disputed 
( Maua s et alJll990l) . the existence of some cool CO in the at- 
mosphere was accepted, though temperature inhomogeneities 
were suggested to be due to the effects of localised mechanical 
heating rather than CO radiative cooling. The presence of cool 
CO in the chromosphere was put beyond doubt four years later 
with the discovery of off-limb CO emission indicating temper- 
atures of less than 4000K in som e places dSolan ki ffi al,|l994). 

Using imaging spectroscopy, Uite nbroek et alJ dl994) stud- 
ied spatial and temporal variations in CO absorption at disk 
centre, observing temperature inhomogeneities of up to 600K 
correlated with hydrodynamic perturbations. This result sup- 
ported the existence of the temperature bifurcation, though 
driven more by dynamic processes than CO cooling and 
magnetic heating as suggested originally. This was picked 
up on in a subseque nt revision of the bifurcation model 
( lAvres & Rabin 1996J) in which the cool component was 

du bbed the 'COmosphere' . 

ICarlsson & Steinl d 19951 Il997l) demonstrated that chromo- 
spheric Can H2V and K2V grains are formed by acoustic 
shocks, using a ID non-LTE hydrodynamic model. Their re- 
sults s trongly supported the notion pushed by Avre s & Rabinl 
d 19961) of a bifurcated atmosphere with differences caused by 
dynamic effects, even also suggesting some influence by CO 
cooling owing to the inclusion of instant aneous chemical equi- 
libri um (ICE) CO o pacities in the model. Kalko fen et al.l dl999) 
and iKalkofenl d200 ll) argued against the Carlsson-Stein model 
because it did not adequately represent the entire solar en- 
ergy output, and so could not reproduce the observations of 
ubiquito us chrom ospheric emission. Howeve r, this was ne ver 
ICarlsson & Steinl s inte ntion, as discus sed bv lAvres| d2002l) in 
an extensive rebuttal of Kalk ofen et all s criticism. |Ayres went 
on to show that despite not being designed with CO diag- 
nosti cs in min d, the Carlsson-Stein mode l far outp erformed 
Kal kofen et al.l s preferred model (VAL; IVernazza et aljfl973l 
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19811) wh en put to work on sol ar CO lines. 
Meanwhile, Uiten hroekl d200 0b) performed detailed ID 
NLTE simulations of CO absorption in the solar disk, re- 
confirmi ng that the LTE app roximation is appropriate in this 
context. llJitenbroekl d2000al) also studied 3D LTE CO line 
formation using a singl e snapshot from an early version 
dStein & Nordlundl 1 19891) of the current state-of-the-art 3D 
atmosp heric models dStein & Nordlundl Tl998t lAsplund et alJ 
2000]J. The calculated line profiles were deeper than those 



observed by the ATMOS mission, something llJitenbroekl at- 
tributed to the ICE ap proximation. 

Growing concern dlJitenbroekll2000albl: lAvresl f2002) over 
the appropriatenes s of the ICE approximation was addressed 
by lAsensio Ramos et alJ d2003l) . who investigated departures 
from ICE in CO lines using the Carlsson-Stein model. Results 
showed significant concentrations of cool CO gas extending 
just 700 km above the Sun's surface, with the ICE assumption 
valid belo w this height. 2D s imulations with non-equilibrium 
chemistry ( Wedemever-Boh m et alJl20 05) have confirmed that 
the ICE approximation is valid in the photosphere and low 
chromosphere. Despite employing non-equilibrium chemistry 
in both the hydrodynamic and radiative transfer simulations, 
modelled CO lines were still considerably deeper than ob- 
served profiles. This is probably because the temperature struc- 
ture of the excerpt from their 2D model is too cool, producing 
an effective surface temperature more than 300K lower than 
the true solar value. Furthermore, 2D simulations tend to give 
erroneous photospheric line profiles compared with 3D calcu- 
latio ns because of t he geometric restriction on the convective 
flow dAsplund et all2000al) . 

1.3. Solar C and O Isotopic Abundances 

The solar carbon a bundance has recently been revised 
dAsplund et alJl2005al) fro m the log e c = 8.52 ± 0.06 given by 
ICirevesse & Sauvall dl998l) Analyses of ICil Ci CH and C 2 
dAllende Prieto et alJ I2OO2I lAsplund et ail l2005bl) now indi- 
cate that the logarithmic solar carbon abundance is in fact 
8.39 ± 0.05. These results provide a comparison for the solar 
carbon abundance determination we carry out using CO lines. 
Considering the difficulty of reproducing observed CO lines 
with theoretical calculations, these provide a very stringent test 
of the model atmosphere. Asymmetries of the CO lines con- 
stitu te an even great er challenge for theoretical atmospheres 
dBlomme et alll994l) . 

Early measurements of the solar 12 C/ 13 C ratio using CH 
lines were typ ically higher than t he terrestrial reference ratio 
of 89.4 ± 0.2 (Coplen et al. 2002), and exhibited large uncer- 
tainties (see lHarris et al.ll 19871 for an inventory of early mea- 
sures). The most reliable (lowest uncertainty) determinations 
to date utilised CO lines, produ cing near-terrestri al ratios of 
84 + 9 dHallll 19731) and 84 ± 5 dHarris et al.lll987l) . The solar 
16 0/ 18 ratio was also measured, producing values of > 500 
and 440 ± 50 respectively. These are simil ar to the re presenta- 
tive terrestrial ratio of 498.7 + 0.1 dCoplen et al.120021) . 1 



1 The geological manner of expressing such measurements is in 
'permil', a deviation in parts per thousand from the terrestrial value: 
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where X is the element in question, a refers to its reference isotope 
and b to the isotope in question. It should be realised that the "terres- 
trial ratios" given above are from standard substances chosen to define 
the zero-point of this scale; in reality, the terrestrial isotopic composi- 
tion vari es far m ore than t he errors attached to these values appear to 
suggest iRosman & TavloJl998l : fCoplen et all2002l) . 
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Different lunar regolith (surface) analy ses of the solar 
wind have indicated 6 n C < -105 ± 20 ( Hashiz ume et alJ 
|2004 and 30 < £ 1 3 C < -30 (combined results reviewed 
in IWiens et aljEoO^ . 5 18 values tabulated in the same re- 
view for direct measurements of solar wind particles by two 
different satell ites give 110 + ^" and 120^ 28 [J. For comparison, 
the results of Ha rris et alJ 111987ft on the same scale imply 
(5 13 C = 59^ and 6 ls O = 130+}^, without the inherent uncer- 
tainty in the scale taken into account. The most recent lunar re- 
golith analyses differ in their estimates of solar £ 18 Q, with one 
suggestin g a negative value ( Hashizume & Chaussidon 2005, 
see also Davis 2005) and the other producing <S 18 = 50 
( Irel and et alJ l2006). Fractionation that might occur in th e so- 
lar w ind is still not properly quantified however JWiens et alJ 
120041) . so such results may say as much about this process 
as they do about the actual photospheric ratios. Solar chem- 
ical evolution is also thought to effect abundance agreement 
between the current photosphere and lunar inclusions irradi- 
ated by the prehistoric solar wind, though the predicted dif- 
ference over the Sun's lifetime is less than 10 perm il in <5 13 C 
or 5 18 Q jTurcotte & Wimmer-Schweingruberl20 02). Recently, 
lYurimoto & Kuramotol J2004 see also lYinll2004 proposed a 
model of solar system formation whereby oxygen isotopic frac- 
tionation occurs in the protosolar cloud due to UV irradiation. 
This theory predicts a solar £ 18 of -50, something that can be 
directly tested with the 18 0/ 16 measurement to be performed 
herein. 

1.4. This Study 

In this paper, we investigate CO absorption line formation 
in the solar spectrum using a three-dimensional radiation- 
hydrodynamic simulatio n of the Sun ' s atmo sphere. This anal- 
ysis is similar to that of lUitenbroekl J2000al) . but with an im- 
proved 3D model as its basis. This enables a better description 
of CO and a more precise evaluation of the reality of the model 
in the context of CO, as well as simulation of line formation in 
the disputed COmosphere. It should be noted that the model is 
not intended to reproduce the ubiquitous chromospheric emis- 
sion nor near-limb CO lines, as it does not include magnetic 
fields nor extend to great enough heights. 

Using the same model and data, we also arrive at new esti- 
mates of the solar carbon abundance, 12 C/ I3 C and 16 C/ ls O ra- 
tios using appropriate sets of weak 12 C 16 0, 13 C 16 and 12 C 18 
lines. Unfortunately, the 12 C n O lines clearly identified in the 
ATMOS spectra are too weak and perturbed to permit any de- 
fensible estimate of the solar 16 0/ 17 ratio, in our opinion. The 
isoto pic ratios co nstit ute an improvemen t over the ID analy- 
ses of lHalll Jl973h and lHarris etaD Il987b . thanks to the use of 
3D model atmospheres, more recent gf values, more appro- 
priate line lists and better observations. We believe that our 
ratios are also supe rior to alter native determinations carried 
out very recently bv lAvres et ail |2006), as will be detailed in 
Sect.0 The determination of the solar carbon abundance using 
CO lines prov i des an important comparison with the study of 
lAsplund et ail J2005bl) based on [C 1], C 1, CH and C 2 lines. It 
also circumvents past trouble in using CO lines and ID mod- 



els to measure the solar carbon abundance (e.g. lGrevesse et al. 
1 1 9951) . 

The 3D model and spectral line lists analysed are briefly 
described in Sect.|2]and Sect, ^respectively, whilst Sect.[4]de- 
tails our manipulation of the ATMOS data. The study of the 
detailed shapes of strong CO absorption lines and their impli- 
cations for the temperature structure in very high atmospheric 
layers is presented in Sect. [5] In Sect. [6] we present abundance 
and isotopic ratio determinations, based on weak CO lines. 
Comparisons with earlier work are given in Sect. [7J and our 
conclusions are summarised in Sect. [8] A detailed description 
of the apodization procedure applied to modelled spectra is 
given in AppendixIXl AppendixIBlprovides derivations of scal- 
ing factors used in the analysis, and Appendix ICl contains full 
line lists. 

2. Model Atmospheres and Line Formation 
Calculations 

The hydrodynamic simulation used was described by 
Aspl und et alJ J200 0b). Specifically, it covers a physical area 
6 x 6 x 3.8 Mm of which about 1 Mm is above the optical 
solar surface (tr oss w 1), at a resolution of 200 x 200 x 82. 
The simulated domain was bounded below by a transmit- 
ting boundary, and above by an extended transmitting bound- 
ary across which the density gradient was kept hydrostatic. 
Horizontal bound aries were periodic. The MHD equa tion- 
of-state jHummer & Mihalasl Il988l iMihalas et alJll988h and 
Uppsala opacities were used, with LTE assumed. Continuous 
and line opacities were calculated using opacity binning. The 
model included no tuneable free parameters, and was charac- 
terised as solar by the accepted gr avity, effective temperatur e 
and standard solar composition of lGrevesse 
A numerical viscosity was employed to stabilise the simula- 
tion. About 100 snapshots of the convective simulation were 
stored, representing approximately 50 min of real solar time. 
Radiative transfer calculations for spectral line formation were 
carried out over an interpolated 50 x 50 x 82 grid, using the 
Uppsala equation-of-state and opacities with molecular densi- 
ties determined under the ICE approximat ion. LTE was a gain 
assumed, as validified for the CO lines by Ayres & Wiedemann 

j — ~n 1 1 j 1 ^ 

(1989) and Uitenbroek (2000a). This is the same atmosphere 
and line formation code recently used to re-derive the solar 
abundances of all elements between Li and Ca (e .g. lAsp lund 
120001 |2004 lAsplund et alJbOOOd. Eooi 12005 albt Asplund et 
al., in preparation;|Allende P rieto et all200lLl2002l) . 

For comparative purposes, abundance calculations were 
carried out using four different model atmosphere s: the 3D 
hydro dynamic model, the ID HM m odel l|H oIweger & M ullet 
1 1974. the ID m arcs model JGustafsson et alJ 1 19751 
Asplu nd et alJ Il997l> . and a contraction of the 3D model 
into the vertical dimension only, which we designate '1DAV' 2 . 
The horizontal averaging used to produce the 1DAV model 
was performed over surfaces of common optical depth rather 
than geometrical height, as the optical depth scale has most 
relevance to line formation. All three ID models included 

2 'One Dimensional AVerage' 
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Table 1. Oxygen and nitrogen abundances used in line forma- 
tion calculations for each mod el in th is study. 3D, HM and 
marcs values are from Asnlu nd et ail 12 004). ID AV ox ygen 
calculated using the same lines as Asplu nd et alJ (|2004) and 
IDAV nitrogen abundance simply set to the 3D value. 



Model 


3D 


HM 


MARCS 


IDAV 


log 6 


8.66 


8.85 


8.80 


8.70 


loge N 


7.80 


7.95 


7.90 


7.80 



E 

2 5000 



-I -;rq " CO 



hv = 2 "C.0 



"C"0 



microturbulence of lkms -1 . HM and marcs were largely 
insensitive to microturbulence, with abundances reduced by 
only 0.01-0.02 dex going from £ t = to £ t = 1 kms _I , 
whilst the IDAV model exhibited an increase of 0.06-0.07 
dex in abundance with £ t = 1 km s instead of £ t = 0. 
Microturbulence was not fitted for in the ID models, as 
they were used simply for comparative purposes rather than 
highly accurate determinations. As it turned out, the model 
most sensitive to microturbulence (IDAV) exhibited the least 
trends in equivalent width and excitation potential anyway. 
The oxygen and nitrogen abundances used for the 3D, HM 
and marcs models were indi cative of the most rec ent values 
produced with each model jAsplund et al.l Eo05al where in 
the case of oxygen, values are from vibration-rotation and 
pure rotational OH lines only). The oxygen abundance used 
for the IDAV model was derived self-consistently from the 
same OH lines used to set the oxygen abundanc es of the 
other models, and employed previously by Asp lund et alJ 
(2004). The nitrogen abundance for the IDAV model was 
simply approximated as the 3D value in the absence of any 
appropriate line calculations. Oxygen and nitrogen abundances 
used in the line calculations with the different models are 
shown in TableQ] 

3. Spectral Lines 

We utilise six disti nct CO line lists, w ith gf values calculated 
from lGoorvitch & Chackerianl i 1994. We ado pted the 12 C 16 
dissociation energy of lEidelsberg et all ill 9871 11.108 eV). All 
lines were selected to be free of blends. Typical formation 
heights (estimated for the marcs model atmosphere) are shown 
in Fig. H CO line formation in higher atmospheric layers is 
examined in Sect. [5] via bisector analysis of 31 12 C ls O lines 
(Table IC.U . These are all strong lines that are formed high 
in the atmosphere. As the height of their formation and hence 
velocity signatures strongly depend on the temperature struc- 
ture of the atmosphere, these lines can probe the disputed 
'COmosphere' region. 

Abundance calculations in Sect. [6] draw upon the five re- 
maining lists, with equivalent widths measured in the ATMOS 
ATLAS-3 spectrum (cf. Sect. @}. The primary investigation is 
performed with a set of 13 weak 12 C 16 lines (Table IC2l of 
high excitation, as well as sets of 16 13 C ie O (Table IC~3t and 
15 12 C ls O (Table IC~4l lines. These lines all have significantly 
lower formation heights than the strong 12 C ls O lines consid- 
ered for Sect. 0(Fig.^. Being CO features, in absolute terms 
they still form reasonably high in the atmosphere, but were se- 



Fig. 1. Approximate optical depths of core formation for the 
different CO line lists used in Sect. |5] (strong 12 C le O) and 
Sect. |6](all others). The temperatur e structure shown is th at of 
the standard 3D model atmosphere (Asplun d et all20 00b). 



lected as forming low enough for abundance determination us- 
ing the 3D model. 

The further two lists consist of 15 low excitation (LE) 
12 C 16 (Table |C5} and 66 first overtone (Av = 2) 12 C ie O 
(Table IC6l lines. The LE lines are formed across a broad opti- 
cal depth range high in the atmosphere (due to the low temper- 
atures required for significant populations in the lower energy 
levels), and the Av = 2 at around the same height as the three 
primary sets (Fig.Q. These supplementary lists are used to de- 
rive alternative carbon abundances (and hence isotopic ratios) 
to the primary (weak) I2 C 16 list, exploring the abundance per- 
formance of the model over higher layers and overtone bands. 
The three different 12 C ie O line lists, which also have different 
temperature sensitivities, have traditionally produced widely 
varying carbon abundance measures using ID models, so a 
consistent result in 3D would greatly increase confidence in 
the new C abundance and isotopic ratios. 

4. The ATMOS Infrared Solar Spectrum 

The ATMOS instrument is a Fourier transform spectrograph 
(FTS). It was carried on a series of space shuttle missions be- 
tween 1985 and 1994, retrieving pure solar and atmospheric 
occultation spectra between 625 and 4800 crrT 1 . Given the in- 
strumental resolution (0.01 cm -1 . lFarmej ri994) and taking a 
representative wavenumber of 2000 cirT 1 (corresponding to a 
wavelength of 5000 nm) for the section of the spectrum pre- 
dominantly considered in this paper, the instrument's resolving 
power for our purposes is: 



2000 



k 

R = — 

Ak 0.01 



200000. 



The solar disk-centre spectrum from the 1994 ATMOS 
ATLAS-3 mission 3 was normalised, with the continuum level 
defined as the highest intensity point in the extracted sec- 
tion. Given the high signal-to-noise of the data (~1000:1 
around 5000 nm by our measurements), detailed averaging 
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was deemed unnecessary. The Sun's gravitational redshift of 
633 m s _1 was removed from the resulting spectrum. The orig- 
inal ATMOS spect rum was apodized using a medium Beer- 
Norton function jFarmegll99 4: Gunson 2004). We hence ap- 
plied the same apodization to our model output; this essential 
procedure is described in detail in AppendixlAl 

To compute CO line bisectors, the absorption profiles were 
split at their centres, with each side separately interpolated to 
a 0.01 normalised intensity unit resolution scale using cubic 
splines. Bisectors were then calculated as the series of wave- 
length midpoints between the two profile halves. Wavelength 
errors were estimated conservatively for each bisector point 
from the previously published (wavelength-dependent) signal- 
to-noise of the data ( Abrams et al. 1996), rather than the higher 
values we measured from the atlas . Bisectors were truncated 
above 0.98 normalised intensity, as closer to the continuum 
they are often dominated by observational noise rather than line 
asymmetry. 

5. CO Line Shapes and Asymmetries 

5.1. Analysis and Results 

In this section, we restrict the discussion to the strongest CO 
lines, formed in the very high atmo spheric laye rs embroiled in 
the current temperature debate (e.g. Avres 2002, and references 
therein). Given their extreme formation heights, such lines are 
not ideal abundance determinants, as our study below demon- 
strates. In Sect. |6j we instead present C and O isotopic abun- 
dances based on weaker CO lines formed in slightly lower lay- 
ers, which our 3D hydrodynamic simulation reproduces more 
faithfully. 

Modelled profiles were fitted in wavelength shift and abun- 
dance by minimising the x 2 likelihood estimator. The wave- 
length fitting was permitted to allow for the possibility of sys- 
tematic errors in the ATMOS wavelength calibration and/or the 
laboratory wavelengths of the lines considered. Unfortunately, 
unfitted modelled line strengths were consistently weaker than 
those observed, requiring that we fit an unrealistically high car- 
bon abundance to produce the best profile agreement. Whilst 
resultant strong CO line profiles showed reasonable agreement 
with ATMOS profiles, their bisectors were of a \-shape (Fig.|2] 
dashed lines) rather than the c-shape observed (solid lines with 
error bars). 

The hydrodynamic simulations have undergone some slight 
improvements since the original model atmosphere was gener- 
ated in 1999, so in an attempt to improve agreement we re- 
generated the CO line profiles using the most modern version 
of the simulations available. The newer code contained an im- 
proved treatment of the numerical viscosity, and drew on an 
updated Uppsala package containing slightly improved con- 
tinuous opacities. The version of the MHD equation-of-state 
was also more recent, including argon in addition to the 16 
other most abundant elements, as well as drawing on slightly 
altered abundances. The manner in which radiation pressure 
and energy were included by the MHD tabulation program was 
also slightly altered. Finally, the mean density and tempera- 
ture structures used in each case to generate equation-of-state 



and opacity tables were slightly different, though it is not clear 
whether this was a primary difference or a secondary effect 
caused by other updates. The final temperature structures were 
very nearly identical, however. 

Inspection suggested that the failure to properly reproduce 
strong lines might have been due to incoming gas from the up- 
permost layer not having time to properly cool and equilibrate 
to its surroundings before it reached CO line formation heights. 
To test this hypothesis, a version of the new atmosphere extend- 
ing to heights of almost 1 .2 Mm was also created. The extended 
simulations were run at a resolution of 50 x 50 x 88 and those 
with the original extension at 50 x 50 x 82. The lower reso- 
lutions were chosen due to computational time constraints, as 
the purpose was not to produce the most accurate description 
of the solar atmosphere available, but to provide a qualitative 
idea of any improvement in CO bisectors due to the later code 
or boundary extension. This is an appropriate approach for the 
qualitative analyses carried out in this Section, but not for e.g. 
abundance analyses, as carried out in Sect. [6] 

Profiles were simulated for all 31 lines, apodized and fit- 
ted. Using the regenerated 3D model atmosphere, very good 
agreement was obtained with observed profiles and bisectors. 
Unfitted line strengths were still smaller than observed, but the 
discrepancy was much less than seen previously. Fitted abun- 
dances were hence still unrealistically high, though the dis- 
crepancy with abundances from weak lines (derived in Sect.|6j 
was only about half that produced by the original model atmo- 
sphere. The scatter in abundances from strong lines was also 
halved when the regenerated model was employed. Some ex- 
ample resultant profiles and bisectors from the (non-extended) 
regenerated model atmosphere are shown in comparison to 
their ATMOS counterparts in Fig.|2] Bisectors derived with the 
older version of the model atmosphere are also given for com- 
parison. The average reduced x 2 value for the agreement be- 
tween observed and theoretical profiles was an order of magni- 
tude smaller using the later model atmosphere. Profiles and bi- 
sectors derived with the extended simulation showed no differ- 
ences to those produced with the corresponding non-extended 
version (data not shown). We hence focus mainly on the results 
from the unextended, later version of the model atmosphere for 
the remainder of our analysis of line shapes and asymmetries. 

5.2. Discussion 

To our knowledge, the only pr ior study of observed asymme- 
tries in solar CO lines is that of Blom meet alJ Jl994l) . Without 
even considering modelled profiles and bisectors, it is plain 
that the observed ATMOS bisectors do not agree with the \- 
shaped bisectors found previously, a significant result in it- 
self. The reason for this discrepancy appears to simply be that 
iBlomme et alJ i 19941) extended their bisectors only to intensi- 
ties of about 0.94 of the continuum intensity, whereas those 
of Fig. El extend t o intensities of 0.98. Given that the bisectors 
of Blom me"et alJ i 19941) were taken from the more noisy SL- 
3 flight of the ATMOS instrument rather than the ATLAS-3 
flight, this lesser extension is not surprising. 



6 



P. Scott et al.: Solar line formation: VII. CO lines and isotopic ratios 




Fig. 2. Example spatially and temporally averaged, disk-centre synthesised strong CO line profiles and bisectors (solid lines 
without error bars), shown in comparison to ATMOS profiles (diamonds) and bisectors (solid lines with error bars). The profile 
agreement is quite good, although as discussed in the text this requires an unreasonably high C and/or O abundance; the bisectors 
based on the new 3D model (see text) show better agreement with observation than those of the original 3D model version (dashed 
lines). The solar gravitational redshift was removed from the ATMOS spectrum, and synthesised profiles have been apodized and 
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The hydrodynamic simulations appear to have been suffi- 
ciently improved for some of the deficiencies in the original to 
be removed. Speculating exactly which minute change or com- 
bination of changes described in Sect. 15.11 is responsible for 
such a higher-order effect as the turning of one end of a bi- 
sector is exceedingly difficult. All that can confidently be said 
on this matter is that it is very unlikely the introduction of ar- 
gon to the equation-of-state had any effect. Beyond this, any 
combination of the tiny changes in opacity and abundance in- 
puts, viscosity and radiation pressure treatments or initial mean 
structure could be responsible. 

However, whether the improvement is attributable to res- 
olutional ef fects rather tha n actual model differences must be 
considered. lAsolund etai1 f2000a) showed that there is an in- 
crease in strong Fe line width going from the resolution of the 
later model to that of the original, with a preferential filling in 
of the bluewards wing of profiles over the redwards wing. This 
is due to the asymmetric way in which low resolution models 
truncate the vertical velocity distribution, removing more of the 
high upwards than downwards velocities. This correspondingly 
manifests itself in a bluewards turning of bisector tops with 
higher resolution, or equivalently, as a redwards 'turning-back' 
of bisectors with insufficient resolution. 

This 'turning-back' has certainly occurred in the current 
study, so the important question is not whether, but to what ex- 
tent the improved bisectors and profil es are attributable to the 
lower resolution. The bisectors o f lAsolund et ail (|2000a ) show 
a maximum redshifting of the bisector top of 1 50 m s~ 1 with the 
change from a 200 x 200 x 82 to 50 x 50 x 82 model. They 
also show a redwards shifting of bisector feet of approximately 
100 m s _1 , resulting in a net redwards movement of 50 m s _I of 
the bisector top relative to the foot. In contrast, the newer model 
bisectors demonstrate a bare minimum of 250 m s _1 net shift of 
the top relative to the foot, compared with the (dashed) origi- 
nals. It would also seem a rather unlikely coincidence for the 
tenfold difference in x 1 values to be due purely to resolutional 
effects. The improvement of profile shapes with the regenerated 
model over the original version is also borne out by a much im- 
proved convergence in fitted line depths and strengths. Fitted 
abundances from the 3 1 lines show far less scatter in the case of 
the later model, and actual abundance values are also closer to 
those derived later using weak lines. These changes support the 
notion that resolutional effects do not play a dominant role in 
the bisector and profile improvements, as the effects of reduced 
resolution are expected to b e the opposite of that observed, i.e. 
an increase in scatter ( Asplund et al. 2000a) and abundances 
derived from CO lines (this paper, Sect. 16. 4> . On the basis of 
these differences, we estimate the contribution of resolutional 
effects to be approximately half of the improvement seen. 

Whilst agreement has been improved, the final results indi- 
cate that the description of the uppermost layers in the model 
atmosphere is still not perfect. This is not surprising, as many 
of the physical approximations used begin to break down in 
the highest layers. Processes such as non-equilibrium hydro- 
gen ionisation, non-LTE continuum formation and extended 
CO cooling related to non-equilibrium chemistry might actu- 
ally be required for an accurate description of such regions. 
Ignoring these non-equilibrium processes could lead to slightly 



anomalous velocity or temperature structures in the top layers, 
issues discussed in more depth in Sect. 16.41 Indeed, it turns out 
that the most likely problem is a slight overestimation of tem- 
perature in the uppermost layers of the model, possibly more 
so in intergranular lanes than granules. This is indicated by the 
overly high abundances derived from strong and LE lines (as 
will be seen in Sect. 16.41 . as well as our own more extensive 
internal investigations into the bisector discrepancies and the 
regeneration exercise (data not shown). Though imperfect, the 
agreement that was seen with observed bisectors does however 
indicate that the 3D model at least captures some of the essen- 
tial physical features of these very high layers, which is encour- 
aging. 

5.3. Temperature Structure and CO Distribution 

To use a model to draw accurate conclusions about the tem- 
perature structure of the solar atmosphere and the distribution 
of CO within it, one would ideally require that synthesised 
CO lines demonstrate excellent agreement with observations 
from very high atmospheric layers. Unfortunately, despite be- 
ing the best effort to date, this has not been the case in this 
study. However, considering that we were able to at least par- 
tially reproduce such difficult indicators as the strong line bi- 
sectors, and more importantly, that all our investigations point 
towards the mean temperature structure in the upper layers of 
the model atmosphere being slightly too high, we can use the 
existing 3D model atmosphere to indicate temperature upper 
bounds for the disputed layers. The temperature structure of 
Fig. ^ hence suggests that cool gas does indeed exist in the 
lower chromosphere, with an average temperature of <4000K. 
Such an upper bound is consistent wit h the extended 3D LTE 
atmosphe re of | Wedemever et alJ J2004I) . and the 3500K pro- 
posed by Avres (2002). Furthermore, the temperature at the 
site of the previous 'minimum' also seems lower than previ- 
ously thought, at <4200K rather than -4500K. Owing to the 
inhomogeneity of the model utilised here, we also tentatively 
suggest the existence of some persistent gas with T < 3700K at 
COmospheric heights, and even intermittent gas temperatures 
of less than 2000K. 

The lack of any profile or bisector improvement in the ex- 
tended simulations would seem to indicate that the extra lay- 
ers above 1 Mm probably do not contribute to CO line for- 
mation at disc centre for the lines considered here. This also 
suggests that the downflowing gas entering the upper layer 
of the non-extended simulations does not play a role in the 
line formation until it has passed significantly beyond the 
top boundary layer and been permitted to adjust to the tem- 
perature of its surroundings. Hence, gas above a height of 
~0.75 Mm might be identified as also generally not contribut- 
ing to disk-centre solar CO line formation (except perhaps 
for some extremely strong lines). Furthermore, as it will take 
some distance to achieve temperature equilibration once gas 
has passed out of the topmost layer of the simulation, our re- 
sults could be seen to suggest an uppermost effective extent of 
the COmosph ere of around 700 km, c onsistent with the indi- 
cations of lAsensio Ramos et alJ J2003t) . It should be noted that 
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this would not however preclude the existence of cool(er) gas 
at greater heights, just suggests an uppermost extent of CO in 
significant concentrations. 

6. C and O Isotopic Abundance Determinations 

6.1. Dealing with CO and Isotopes 

Working with a known oxygen abundance, a fitted abundance 
difference for CO lines can be approximately interpreted as 
a difference in solar carbon abundance, even though strictly 
speaking it is actually indicative of a difference in CO abun- 
dance at the height of line formation. In order to use CO lines 
to derive a solar carbon abundance, the input carbon abundance 
was therefore iteratively altered according to the average fitted 
abundances of the CO lines in the previous iteration, until con- 
vergence was obtained. Oxygen concentration effects the equi- 
librium position of the CO formation-dissociation reaction, so 
carbon abundances derived in this manner are dependent upon 
the adopted oxygen abundance. Input oxygen abundances were 
hence carefully chosen (cf. Sect. |5J and kept constant through- 
out the study. 

Using the above technique, the solar 12 C abundance was de- 
termined separately using the three different 12 C 16 line lists. 
Since the 3D model input abundances do not differentiate be- 
tween isotopes, these determinations were performed by alter- 
ing the total input carbon abundance, though the final abun- 
dance arrived at is not overall abundance, but the carbon abun- 
dance if all C in CO were contained in l2 C l6 0. This can ap- 
proximately be called the 12 C abundance, even though it does 
not include any 12 C tied up in 12 C n O or 12 C ls O. Though they 
are minimal given the high values of 16 0/ 18 and 16 0/ 17 0, 
these contributions were taken into account in the calculation 
of a total carbon abundance by the use of a fractional scalefac- 
tor (Appendix lB.3> . 

The concentrations of 13 C 16 and 12 C ls O relative to 
12 C ls O are identical to the 12 C/ 13 C and 16 0/ 18 ratios. In order 
to determine the concentrations of the CO isotopomers, given 
the lack of provision for isotopic differentiation in the model, 
the isotopomeric lines were treated in the radiative transfer sim- 
ulations as if they were created by 12 C 16 0. In order to do this 
correctly, the radiative transfer code had to be altered to include 
provision for a mass scalefactor and an opacity scalefactor (re- 
fer to AppendixlBl. 

6.2. Abundance Calculations 

For each iteration, three profiles with log gf values differing by 
0.2 dex were calculated for a given line and interpolated be- 
tween using cubic splines to arrive at the desired line strength. 
As described in Sect. |2] abundance calculations were carried 
out for all five line lists with four different mo dels: 3D, HM, 
marcs a nd 1DAV. The 3D model used was that of As plund et alJ 
(2000tJ, the same model atmos phere used for all previous 
3D abundance determinations (e.g. lAsplund et all2000d.l20o4 
2005a b; Asplun dl2000l Eb04). Since accuracy is paramount in 
abundance determinations, the regenerated version described in 
Sect. |3 would have been inappropriate given its reduced reso- 



lution. This model was also used over the regenerated version 
for consistency and comparability with the previous abundance 
determinations. 

Both equivalent width and profile fitting (the latter via % 2 - 
analysis) were tested upon the first iteration of the weak 12 C ls O 
lines, with virtually no difference found in calculated abun- 
dance. However, some of the lines could not be effectively 
profile fitted by minimising the x 2 statistic without some form 
of masking, due to the presence of other nearby lines in the 
ATMOS spectrum. For this reason, equivalent width fitting was 
used for all subsequent abundance measures, as this problem 
was only likely to worsen when the weaker isotopomeric lines 
were considered. The same Beer-Norton apodization as used in 
Sect. [5] (medium BNA, characteristic velocity 1.5 kms -1 ) was 
also applied to each line, though more for the sake of consis- 
tency than anything else. Convolution with any function nor- 
malised to unit area (as the medium Beer-Norton used was) 
should be an area-preserving operation, thereby not effecting 
line equivalent widths nor therefore the abundances determined 
with them. For completeness however, it should be noted that 
in the case of the Av = 2 lines, a BNA characteristic velocity 
of 0.75 kms~' was used, reflecting the higher resolving power 
of the ATMOS instrument at the shorter wavelengths of these 
lines (see Table lC.ol . 

6.3. Results 

Derived abundances are plotted as a function of equivalent 
width and excitation potential for each line in the weak 12 C ie O 
(Fig. LE 12 C 16 (Fig.EJ, Av = 2 I2 C le O (Fig.0, 13 C ie O 
(Fig ■ EJ and 12 C ls O sets (Fig.0. The only appreciable trend 
evident in the 3D results is with equivalent width in the case of 
the LE lines, though these lines do also produce a smaller trend 
in excitation potential in 3D. The HM and marcs models dis- 
play trends in equivalent width and excitation potential for the 
12 C 18 lines, and excitation potential in the Av = 2 lists. In ad- 
dition, the marcs model exhibits a trend with equivalent width 
in its LE results. In the case of the 1DAV model, no significant 
trends are seen except for a small slope with equivalent width 
in the 12 C 18 list. 

Final abundance measures and isotopic ratios as defined by 
the weak, LE and Av = 2 12 C 16 lines are tabulated (Table|2}. 
Errors in 12 C abundances are given by single standard devi- 
ations within individual 12 C 16 line lists. Errors in bulk C 
abundances are single standard deviations within each 12 C 16 
list when a fractional scalefactor (cf. Appendix IB. 31 was em- 
ployed, as uncertainty in the scalefactor due to isotopic abun- 
dance uncertainties was negligible compared to line-to-line 
scatter. Similarly, errors in isotopic abundances are a single 
standard deviation within the relevant isotopomeric list. Errors 
in isotopic ratios were determined in logarithmic units as the 
sum in quadrature of errors attached to logarithmic 12 C and 
isotopic abundances. Asymmetric errors in absolute ratios are 
given as corresponding deviations due to these calculated (sym- 
metric) logarithmic errors. Overall results in 3D from different 
indicators are summarised in Table[3] 
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Table 2. Abundances and isotopic ratios implied by the weak (Av = 1), LE (Av = 1) and Av = 2 12 C 16 line lists as indicator 
of 12 C abundance. Note the large differences between ID and 3D isotopic ratios. Note also that the HM and marcs ID models 
indicate higher abundances than the 1DAV, which in turn produces a higher abundance than the 3D model. 



12 C 16 Lines 


Model 


3D 


HM 


MARCS 


1DAV 


Terrestrial 


Weak 
(Av = 1) 


loge c 
log ei2 C 

log 613 C 

log as 
12 C/ 13 C 

L6 0/ 18 


8.40 ± 0.01 
8.39 ±0.01 
6.44 ± 0.02 
5.97 ± 0.03 
88.8!" 

490!^ 


8.60 ± 0.01 
8.59 ±0.01 
6.76 ± 0.02 
6.31 ±0.04 
68.6! 32 
344! 3 » 


8.55 ± 0.02 
8.55 ± 0.02 
6.70 ± 0.02 
6.25 ± 0.04 
70.6^ 3 
350!^ 


8.46 ± 0.01 
8.45 ± 0.01 
6.58 ± 0.02 
6.12 ±0.03 

74.o ;: 

376! 3 ^ 


89.4 ± 0.2 
498.7 ±0.1 


(Av = 1) 


loge c 
log ei2 C 
log eo c 
log eis 
12 C/ 13 C 

16 0/ 18 


O A O i A A/1 

8.4a ± 0.04 
8.47 ± 0.04 
6.44 ± 0.02 
5.89 ±0.03 
107.6+j 3 
594^ 


o i n no 
8.00 ± 0.02 

8.66 ± 0.02 
6.76 ± 0.02 
6.24 ± 0.04 

80 - 3 -54 

403!« 


8.0 / ± U.UJ 
8.67 ± 0.03 

6.70 ± 0.02 
6.13 ±0.04 

463 +57 


o /in i n no 

8.49 ± 0.02 
6.58 ± 0.02 
6.09 ± 0.03 
80.0!|> 

4orf 5 


89.4 ± 0.2 
498.7 ±0.1 


Av = 2 


logec 
log ei2 C 
log ei3 C 
logeis 
12 C/ 13 C 

16 0/ 18 


8.37 ±0.01 
8.36 ±0.01 
6.44 ± 0.02 
6.00 ± 0.03 

82.8!" 
457 +39 


8.69 ± 0.02 
8.68 ± 0.02 
6.76 ± 0.02 
6.22 ± 0.04 

S3.9 

42 1 +43 


8.58 ± 0.02 
8.57 ± 0.02 
6.70 ± 0.02 
6.23 ± 0.04 

-5.0 :; 

372!^ 


8.51 ±0.02 
8.50 ± 0.02 
6.58 ± 0.02 
6.07 ± 0.03 

82.8!^ 
42L 3 5 8 


89.4 ± 0.2 
498.7 ±0.1 




5 6 7 

Equivalent Wid;h [pm] 




2 3 4 

Equivalent W'dth [prr] 




2.2 2.4 2.6 

Lower Excitation Potential [eV] 



8.75 
« 8.70 
o 



8.65 
8.60 
, 8.55 
8.50 



8.45 



• HM 

o MARCS 



o o 
- -o- 



O o 



* • 

o 



O O o 



2.0 2.2 2.4 2.6 2, 

Lower Excitation Potential [eV] 



Fig. 3. Solar carbon abundances indicated by the weak 12 C 16 lines, displayed according to equivalent width (top) and excitation 
potential (bottom). On the left, filled circles indicate 3D results and open circles 1DAV results, whilst on the right filled circles 
are HM values and open circles marcs results. Trendlines are produced as linear fits to data sets using a minimised^ 2 method 
placing equal weight on each point, with solid lines corresponding to filled circles and dashed lines to open circles. No significant 
trends can be seen in the output of any model. 
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consistent with the abundances derived using the weak and Av = 2 12 C 16 lines. The agreement deteriorates with increased LE 
12 C ie O line strength and therefore formation height. Less prominent trends are also evident in excitation potential for all models. 



6.4. The Solar C Abundance 

The absence of significant trends in 3D model abundances over 
the primary (weak) and overtone (Av = 2) line lists (Figs. [3] 
and|5Jl, as well as the similarity of the resulting average abun- 
dances (log ec = 8.40, 8.37), give us confidence in the accuracy 
of these values. They also suggest that the 3D model is accu- 
rate in the quite high atmospheric layers in which these lines 
form. The agreement of abundances derived using the normal 
weak and overtone lines is quite an achievement, as these indi- 
cators have not generally produced consistent results. The HM 
model for example produces abundances differing by 0.08 dex 
between the two lists (and by nearly three times this between 
CO lines and other in dicators of carbon abundance stated in 
Asnl und et al.ll200.5bl) . It also exhibits a clear trend with ex- 
citation potential for the Av = 2 lines (bottom-right of Fig.|5j, 
further suggesting the inadequacy of its description of overtone 
CO. The same trends are true of marcs, though it at least pro- 
duces consistent abundances between the two line lists. That 
the 1DAV model abundances in these two cases are much closer 
to the 3D results, though still not quite as low, reflects the fact 
that both the mean temperature structure and the temperature 
inhomogeneities play a role in producing different CO-derived 
C abundances in 3D than with HM or marcs. 

The same success was not seen in LE results, where a strik- 
ing trend with equivalent width is present in the 3D case (top 
left of Fig. [4}- The reason for this is very likely the larger range 



of formation depths of these lines and increasing inadequacies 
of the 3D model atmosphere in higher layers. The weaker of the 
LE lines can be seen to produce abundance measures very close 
to the primary and overtone results, as they form at similar 
heights (cf. Fig.0. However, the derived abundance increases 
with greater equivalent width and line formation height, result- 
ing in an average abundance quite a bit higher than the primary 
or overtone diagnostics. 

A number of reasons could be postulated for this discrep- 
ancy. The failure of the chemical equilibrium approximation 
(ICE) with height in the atmosphere would mean that CO den- 
sity was being overpredicted in the models, producing stronger 
lines than otherwise would result from a particular abundance. 
Hence, less abundance would be required to reproduce a given 
line profile than in the non-ICE case, not more as is seen here. 
The breakdown of LTE at height in the atmosphere could possi- 
bly cause overestimated abundances here, though it seems un- 
likel y given the repeated conclusion that CO lin es form in LTE 
jAvres & WiedemanrJl989HUitenbroekj20 00b). Another pos- 
sibility could be that temperature contrast in the upper layers of 
the model is too low, as an increase in this contrast would pro- 
duce lower temperature cool regions, which contribute more to 
increasing line strength than hotter regions would to decrease 
it, due to the increase in CO and the nonlinear temperature de- 
pendence of line formation. However, this seems unlikely as 
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the 12 C ls O lines. Significant trends can be seen in equivalent 



Table 3. Summary of carbon abundances and isotopic ratios produced with the 3D model by the three different 12 C ie O line lists, 
as well as the adopted values. Note the poor agreement of the LE results with the other two, reflecting the trend seen in Fig.0] 
Adopted values were calculated via 2:0:1 weightings of the weak:LE:Av = 2 lists. 



12 C 16 List 


Weak 


LE 


Av = 2 


Adopted 


Terrestrial 


logec 


8.40 ±0.01 


8.48 ± 0.04 


8.37 ±0.01 


8.39 ± 0.05 




12 C/ 13 C 


88.8!" 


107.6!]^ 




86.8+13 


89.4 ± 0.2 


16 0/ 18q 


4901^ 


594:*» 




470+29 


498.7 ±0.1 



our bisector investigations hinted at too much temperature con- 
trast high in the atmosphere. 

Also possible is that the velocity structure of the upper at- 
mosphere is slightly incorrect in the 3D simulations. The marcs 
model exhibited a similar equivalent width trend in implied 
abundances to the 3D model for the LE lines, also suggest- 
ing problems in its upper layers. However, even the strongest 
LE lines were not sensitive to the microturbulence parameter 
used with the marcs model, making this explanation unlikely. 
Similarly, other controls revealed that the problem was not due 
to the adopted collisional damping parameter. The most likely 
explanation for the pronounced trend seen in the 3D LE results 
is that the mean temperature structure at height is slightly too 
high, causing a reduction in line strength for a given abundance 
and conversely, increased abundances for given line strengths. 
This is consistent with the lack of trend but higher average 
abundance indicated by the 1DAV model, as the horizontal 
averaging would smear out the effects of overly hot localised 
regions (e.g. intergranular lanes) at height but produce higher 



abundances overall due to the complete absence of temperature 
inhomogeneities. This is also in line with the bisector discrep- 
ancies seen using strong CO lines. 

The regenerated model of Sect. [5] was trialled on the LE 
and weak 12 C 16 lines to see if agreement could be improved, 
and also to check that the newer model did not produce signifi- 
cantly different abundances where the older model was thought 
to be accurate (i.e. the weak lines). Derived abundances in both 
cases were approximately 0.03 dex larger, virtually identical 
to the effect of low re solution upon iron abundances found by 
lAsolund etailfeOOOal) . 

Despite the poor performance of the 3D model for the LE 
lines, none of the ID models managed to derive much better 
agreement between these and the weak 12 C 16 lines (which 
are regarded as the best CO indicators of carbon abundance). 
In the sense of the LE lines' disagreement with other diagnos- 
tics, the 3D model simply presented no improvement, rather 
than any loss in performance over HM and marcs. In all three 
diagnostics, the carbon abundance in 3D is considerably lower 
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than indicated by HM o r marcs, as has generally been fou nd for 
other species (e.g. lAsplund et alJl2000ct l2004 12005 albh . This 
is a general consequence of extending atmospheric simulations 
to three dimensions, due to the permission of temperature inho- 
mogeneities and the nonlinear temperature dependence of line 
formation. That all three sets of 12 C lfi O lines continued to dis- 
play this effect is a positive comment on the accuracy and con- 
sistency of the 3D model, esp ecially given the past difficulties 
with ID analyses of CO lines jGrevesse et alll995l) . 

Given the relative performance of the different 12 C le O 
line lists and their previously recognised suitability for carbon 
abundance determination, for the final analysis weak lines were 
given a weighting of 2, Av = 2 lines a weighting of 1 and LE 
lines no weighting at all. We note however that abundances 
from the weakest LE lines are in good agreement with those 
from the two preferred line types. Hence, on the basis of the 
lines measured in this study, the final carbon abundance arrived 
at is 

loge c = 8.39 ±0.05. 

The stated error is designed to encompass unqualified system- 
atic errors from e.g. atomic data, the model atmosphere and 
equivalent width measurements; these systematics almost cer- 
tainly outweigh any statistical variability amongst abundances 
indicated by individual lines or lists. 

6.5. The Solar 12 C/ 13 C and 16 0/ 18 Ratios 

The abundances produced by the individual I3 C 16 Oand 12 C ls O 
lines in 3D (left panels in Figs. [6] and [7} show no significant 
equivalent width or excitation potential dependence, hence in- 
dicating no major inadequacies in the relevant atmospheric lay- 
ers of the model. This is also the case for ID calculations with 
the 13 C lfi O lines (right panels in Fig.|6j but not the 12 C lg O lines, 
where trends with equivalent width and excitation potential are 
evident in the HM and marcs results (right panels in Fig. 0. 
The two trends are explicitly linked, as the higher excitation 
lines form lower and therefore have smaller equivalent widths, 
so it is no surprise that one trend exists in equivalent width and 
the opposite is seen in excitation potential. These trends would 
seem to indicate the superiority of the 3D model in the descrip- 
tion of the very weak 12 C ls O lines. 

The derived isotopic ratios in Table [2] using the different 
model atmospheres show a very strong separation between 3D 
and ID results. This is surprising at first, though suggests that 
in this context, the move to three dimensions is an important 
improvement. This is highlighted in particular by the general 
difference between 3D and 1DAV results. Our investigations 
have revealed that, due their lower excitation potentials, the iso- 
topomeric lines are actually over 20% more temperature sensi- 
tive than the weak 12 C ls O lines. At least in the case of the weak 
12 C 16 lines, this is almost certainly the reason for the greater 
decrease in isotopic abundances than the carbon abundance in 
3D. This is because the heterogeneic temperature structure of 
the 3D simulations (which is responsible for reduced abun- 
dances normally) has a greater decreasing effect upon the iso- 
topomeric lines and therefore 13 C and ls O abundances than it 
does upon the derived carbon abundance. Because of this dif- 



ferential temperature sensitivity, the horizontal averaging im- 
plicit in ID models means that they have difficulty reproducing 
the correct solar isotopic ratios. Given their excitation poten- 
tials, the temperature sensitivities of the LE and Av = 2 lines 
are rather similar to the isotopomeric lines. This is the reason 
for the better agreement between ID and 3D isotopic ratios 
from the Av = 2 lines than from the weak 12 C 16 lines. Any 
similar effect in ratios derived from LE lines is lost in the dis- 
tortion of the 3D carbon abundance by the observed trend with 
formation height. 

Because each 12 C ls O list produced different abundances 
and each was compared to the same 13 C ls O and 12 C ls O re- 
sults, different 3D isotopic ratios (and absolute 18 abun- 
dances, since eo was fixed) were obtained. However, the LE 
ratios can be discarded due to the already established invalidity 
of the LE carbon abundance. In the derivation of the final iso- 
topic ratios, the Av = 2 ratios were given a lower weighting in 
the same manner as in the calculation of the final carbon abun- 
dance (i.e. half the weighting afforded the weak ratios). Hence, 
the adopted isotopic ratios are: 



12 C/ 13 C : 



86.8 



+3.9 
-3.7 



16 0/ 18 = 479!| 

Errors are combinations of statistical errors in the weak and 
Av = 2 results. We assume that the systematic errors discussed 
as relevant to the carbon abundance in Sect. l6.4l effect all lines 
considered for ratios equally, so need not be included in final 
ratio errors. This assumption is based on fact that the isotopic, 
weak and Av = 2 lines all form at approximately the same pho- 
tospheric heights, reflecting essentially the same line formation 
processes. 



7. Comparisons with Previous Work 

This study represents a significant step forward fr om the only 
other 3 D investigation of CO line formation to date dUitenbroekl 
2000a). Firstly, we used a more modern 3D hydrodynamic 
model atmosphere code, whereas the previo us study was based 
upon the very early Stein & Nordlund ( 1989) version, running 
at a resolution of just 63 x 63 x 63. Secondly, the current study 
produced profiles spatially averaged over the entire simulation 
domain and temporally averaged over about 100 snapshots cor- 
re sponding to approxim ately 50 min of solar time, whereas that 
of llJitenbroekl (EoOOal) used a vertical slice through a single 
snapshot. Thirdly, the maximum vertical extent of the current 
study was about 1 .2 Mm, compared with 0.6 Mm in the ear- 
lier work. Finally, agreement between the modelled spectrum 
and observation was excellent overall, with reasonable agree- 
ment even found at the level of individual line bisectors, a 
measure highly sensitive to smal l deviations of the simulation 
from reality. As suggested bv lAsensio Ramos et al.l J2003). the 
overly cool version of the 3D model atmosphere employed by 
llJitenbroekl is probably somewhat responsible for the overly 
deep CO line cores he foun d. The recently reduced sola r car- 
bon and oxygen abundances ( Asnlund et al.l2005bl Eo04) have 
also played a role in the success of the cu rrent study, permit- 
ting the deviation of line profiles from the iGrevesse & Sauvall 
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( 1998) abundances in order to improve agreement with obser- 
vation. This option was not realistically available to Uiten broekl 
(2000a) given prevailing wisdom at the time, and in light of the 
current results this probably also contributed to his modelled 
overprediction of CO line depths. 

In the context of spectral line formation, our results also 
constitute an improvement ove r the more recent 2D work of 
Wedemever-Bohm et a This is due not only to the 

better agreement achieved with observation, but also the su- 
periority of our model atmosphere in this context. Apart from 
being 3D rather than 2D and therefore inherently more real- 
istic, the model we employ includes line blanketing, rather 
than treating radiative transfer in the hydrodynamic code as 
grey. In addition, our model exhibits an effe ctive tempera- 
ture in excellent agreement with the Sun's (se elAsplun d et alJ 
2000b), whereas the snapshots Wedemever-Bohm et al. em- 
ploy h ave an effective temperatu re approximately 380K lower. 
Whilst Wedemever-Bo hm et al.l utilise non-equilibrium chem- 
istry, they find that it is not required at the heights of formation 
considered in this paper, justifying our use of the ICE approxi- 
mation throughout. 

The adopted bulk carbon abundan ce is in excellent agree - 
ment with the log e c = 8.39 ± 0.05 of lAsplund etaD d2005bl) . 
which was based on very different indicators (C i, [C i], C2 and 
CH). The current figure therefore firms our belief in the ac- 
curacy of both carbon results and the 3D model atmosphere. 
In comparison, CO results with the preferred lines using the 
HM and marcs models do not agr ee at all with their C 1, [Ci], 
C2 and CH counterparts given by Asplun det alJ J2005bl) : CO 
lines give log e c = 8.60, 8.69 but other lines 8.39-8.53 with 
HM, whilst with marcs CO lines indicate logec = 8.55, 8.58 
and other lines give 8.35-8.46. Given the height of forma- 
tion of the CO lines, this result is a remarkable success for 
the 3D model. The new carbon abundance constit utes a ma- 
jor red uction from the commonly adopted iGrevesse & Sauvall 
( 1998) value of log ec = 8.52 ± 0.06, not t o mention a multi- 
tude of even higher earli er estimates (e.g. lHarris et alJ ll987: 
lAnders & GrevessJ 19891) . 

Using the HM model. lHarris et all i 19871) found isotopic ra- 
tios of 12 C/ 13 C = 84 + 5 and 16 0/ 18 = 440 ± 50, broadly con- 
sistent with the ratios we obtain according to the stated errors. 
In comparison, our own HM results (weighted amongst the lists 
in the same manner as the 3D ratios) were 12 C/ 13 C = 73.4*2? 
and 16 0/ 18 = 368+g. The difference in HM results between 
the two studies evidently reflects some combination of the im- 
proved observations, line lists and molecular data used in the 
present work. 4 However, the effect of these improvements is 
to reduce the ratios, whereas the adopted ratios calculated in 
3D are slightly greater than those previousl y found. Hence, the 
agreement between our results and those of Harr is et a ll lll987l) 
is somewhat fortuitous; any difference between old and new 



adopted ratios primarily reflects the shift from ID to 3D mod- 
elling, with the effects of this shift in fact mostly nullified by 
improved input data. 

We also believe that our ratios are more r eliable than 
those very recently produced independently by lAvres et all 
fcOOfit 12 C/ 13 C = 80 + 3 and 16 0/ I8 = 440 + 20, where errors 
are standar d deviations rather than the lower standard error 
lAvres et all favour - ). They utilise our 1DAV model and a mul- 
ticomponent version of it meant to approximate the true 3D 
case. Discrepancies between isotopic ratios derived in the two 
studies using the 1DAV model can be explained by differ- 
ences in measured equivalent widths and our use of the opac- 
ity scalefactors detailed in Sect. IB.21 5 . We believ e our equiv- 
alent w idth measurements to be superior: whilst lAvres et all 
(2006) utilise Gaussian fits to ATMOS data, we performed di- 
rect integration upon line profiles, since the ATMOS ATLAS-3 
data is of sufficiently high signal-to-noise to make this possi- 
ble and observed profiles always exhib it more exten ded wings 
than a pure Gaussian. Furthermore, all lAvres et alt s Gaussian 
fits to isotopomeric lines are constrained to have a FWHM of 
4.05 kms _I , whereas our integrati on allow ed this value to be di- 
rectly determined. However, Avr es et alJ choose to discard the 
results of our model in favour of an empirical ID model based 
upon a 'Double-Dip' photospheric structure, which has no cor- 
roborating physical or theoretical basis beyond its ability to 
match the observed centre-to-limb variation of CO lines and the 
solar continuum. The derived abundances show marked trends 
in excitation potential and a disagreement between fundamen- 
tal and overtone b ands, which our m odel does not (either in 
this paper or that o flAvres et all2006l) . To compensate for these 
trends, the authors introduce a highly-questiona ble, excitatio n 
potential-dependent change in given gf values. lAvres et all s 
choice of model is made upon the basis of the 1DAV model 
failing to properly reproduce the observed continuum and CO 
core centre-to-limb behaviour. However, the stated discrepancy 
in the continuum centre-to-limb results is about double that 
produced by our own calculations. We intend to return to this 
important issue in a subsequent study, where the temperature 
structure and full results of a large number of observational 
tests will be described. 

The 12 C/ 13 C ra tio we find is in excellent agreement with the 
Coo len et alJ (T2002) terrestrial value of 89.4+0.2. Our ratio cor- 
responds to a 6 l3 C value of 30^ (without taking into account 
the uncertainty in the <5-scale), in agreeme nt with some of the 
lunar regolith measures discussed bv lWiens et al.1 d2004l) . Our 
ratio does not agree wi th the <!> 13 C < -105 + 20 measured by 
Hashizume et al. (12004 . suggesting that theories for the forma- 
tion of organics in the early solar system based upon this low 
value may require revision. We confirm a hig h solar 12 C/ 13 C 
ratio relative to the local ISM (62 + 4: La nger & PenziaJ 19931 



4 Th e log gf values were typically 0.02 dex larger in the lHarris et alJ 
( 1987) study, so these are unlikely to have affected rat ios. None of 
the weak 12 C 16 lines we have used were employed by Harris et all 
Equivalent widths for the few common isotopomeric lines were larger 
on average in their study, leading to an underprediction of isotopic 
abundances and therefore an overprediction of ratios. 



5 Adopting Avres et alT s smaller equivalent widths, and consider- 
ing the resultant average change in abundances over the small number 
of lines common to the two studies, as well as discarding the opacity 
scalefactors, we would obtain 12 C/ 13 C = 80.9 and 16 0/ ls O = 452 for 
the 1DAV model using the weak I2 C I6 list. This is in comparison to 
12 C/ 13 C = 80.7 ± 3.8 and 16 0/ 18 = 443 ± 31 obtained bv lAvres et all 
using the same model. 
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68 ± 15: Mil am et alJ Eool . This may support the notion 
of a galactic enric hment of 13 C over the Sun's lifetime (e.g. 
iMilam etal.l Eo05). Alternatively, the difference could be due 
to greater isotopic inhomogeneity in the galaxy than previously 
thought, and migration of the solar system through it. 

Our derived 16 0/ 18 abun dance is also consist ent with the 
terrestrial ratio of 498.7 ± 0.1 (Coole n et alJl2002h . The corre- 
sponding (5 18 value is 41*5Z (with error in the (5-scale not ac- 
counted for). Despite the lack of error bars on the £ 18 Q _ -50 
prediction of the self-shielding model jYurimoto & Kuram oto 
2004| I Yinl 120041) . it seems unlikely that our results would be 
consistent with this prediction. However, our result is in ex- 
cellent agreement w ith the (5 18 = 50 recently determined by 
llreland et al.l J2006I) from lunar grains irradiated by the solar 
wind. However, given uncertainty about how well the solar 
wind actually reflects the solar composition, and the size of 
our estimated errors, this amazing agreement could be slightly 
fortuitous. 

8. Conclusions 

This is the most successful modelling of CO line forma- 
tion in the solar atmosphere to date. Solar abundance deter- 
minations using CO lines produced a carbon abundance of 
logec = 8.39 + 0.05, and isotopic ratios of 12 C/ 13 C = 86.8^ 
and 16 0/ 18 = 479^. These results represent a significant 
improvement over those of the past due to the combination 
of a state-of-the-art convective 3D model atmosphere, up- 
dated atomic data, better line lists and more accurate obser- 
vations. The carbon abundance is in excellent agreement with 
the recent findings ba sed upon entirely different indicators of 
Asnl und et al.l J2005bl) . suggesting that the past problems with 
CO-derived abundances have been solved. Both ratios are in 
excellent agreement with corresponding telluric values, though 
the oxygen ratio is in even closer agreement with a non- 
terrestrial value inf erred by the latest lunar regolith analysis 
( llreland et all2006T) . Our results also support the existence of a 
COmosphere, with a representative temperature below 4000K. 
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Appendix A: ATMOS Apodization 

Knowing the instrument's resolving power, it becomes possible 
to simulate the observational effect of an FTS. This is simply 
done by convolving modelled spectra with an instrumental line 
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Fig.A.l. The effects of medium BNA upon the same 
4752.2 nm CO line as shown in Fig. [21 The dotted profile has 
had no instrumental profile applied to it, whereas the solid 
line illustrates the effects of medium BNA characterised by 
Act =1.5 kms -1 . Convolution with a simple sine function of 
Act =1.5 kms 1 leaves the line profile virtually identical to its 
unconvolved counterpart. In the interests of clarity this curve 
is not shown. Note the far better agreement with the ATMOS 
spectrum of the apodized than the unapodized profile (reduced 
X 2 measures indicate an order of magnitude difference), vividly 
demonstrating the importance of correctly emulating instru- 
mental effects when working with this data. 

shape (ILS) characterised by the spectral resolution in the ap- 
propriate scale. 

Unfortunately however, FTS instruments are not always so 
well behaved as to exhibit perfect top-hat interferogram pro- 
files. This was certainly the case with the ATMOS ins trument, 
which had an ILS that was far from a perfect sine JGunsonl 
2004). To deal with the uncertainty in the true instrumental pro- 
files of transform spectrometers, apodization is sometimes em- 
ployed. The process of apodization involves the multiplication 
of the interferogram with a profile other than the typical top-hat 
function, such as a Gaussian or simple triangle function, with 
somewhat tapered edges. This of course has the same effect as 
convolving the spectral output with an ILS that is the Fourier 
transform of the apodizing function which one multiplies the 
interferogram by. This smears out the variability in the instru- 
mental profile, and dampens harmonic edge effects produced 
in the resultant spectrum by the finite maximum optical path 
difference of the FTS. 

T he apodiza ti on ap plied by the ATMOS team jFarmerl 
1994; Gunson 2004) was the medium function of 



Norton & Beer I 




, sine a - cos a 



ILS(a-) = 0.26sinca -0.464514- 

a- 

-13.422570 (1 - 3/a2)sinCa 7 +(3/a2)C ° Sa , (Al) 

a 2 

where a — cr is the spectral variable and Act its resolu- 
tion. This fun ction is one of three purposefully created by 
iNorton & BeeJ for their minimal resolutional broadening of 
spectral points and maximal damping of secondary maxima. 
In order to correctly simulate the effects of the ATMOS in- 
strument and the post-processing performed by its operators, 
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we convolved modelled profiles with the same medium Beer- 
Norton function. Model spectra were expressed in terms of 
Doppler velocity, so the convolving function argument was ve- 
locity resolution, i.e. | = 1.5 kms -1 . The striking effect of this 
post-processing is seen in Fig. IA.ll where a predicted 3D line 
profile that has undergone Beer-Norton apodization (BNA) fits 
ATMOS data far better than either its unconvolved counterpart, 
or a profile convolved with a plain sine function. 

Appendix B: Scalefactors 

B.1. Mass scalefactor 

The mass scalefactor allowed the correct molecular mass to be 
used in the line formation calculations, and was simply the ratio 
of the mass of the isotopomer in question to that of the most 
common isotopomer (i.e. 12 C lfl O). These factors, 1.03583 for 
13 C 16 and 1.07 157 for 12 C 18 Q, were ca lculated from nuclear 
weights given bv lAudi & Waostral i 19951) . 

B.2. Opacity Scalefactor 

The opacity scalefactor contained the information about dif- 
ferent isotopomer concentrations and therefore the atomic iso- 
topic ratios, as the line opacities for any particular transition 
in two different isotopomers are proportional to their densities. 
The use of an opacity scalefactor was preferable to simply al- 
tering the input carbon abundance in order to emulate the dif- 
ference in isotopic abundances. This is because the overall car- 
bon abundance was not altered and therefore did not feed back 
on other parts of the simulation like CN line formation or the 
temperature structure, which might in turn indirectly effect CO 
line formation. The opacity scalefactor was necessary anyway, 
in order to include an actual molecular difference (rather than 
just that due to abundances) in opacity between isotopomers. 
This opacity correction arises as follows: 

First, recall that Doppler (thermal) line broadening by a 
species of mass m has a narrow Gaussian characteristic profile 
of the form 



6d(A^) 



1 



,-(A,l/Ai D ) 2 



(B.l) 
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where AA is the distance from the line centre of wavelength Ao 
and 



A 2kT , 
c V m 



(B.2) 



is the Doppler width of the Gaussian. The microturbulent ve- 
locity £ t is an extra 'fudge factor' used in analyses based on ID 
model atmospheres, introduced to emulate Doppler broadening 
due to inherently three-dimensional microturbulence in the gas. 

Now, consider two isotopomers of a diatomic molecule of 
elements X and Y: W = a X c Y and = b X c Y. The opac- 
ity of a species in a given transition depends on the density of 
absorbers, the transition probability and the relevant line broad- 
ening effects. In the case of a weak lines, as those considered in 
this study are, Doppler broadening given by the profile 0d from 
Ea. lB.ll can be approximated as the only significant contributor 
to the latter, meaning that for our molecule S, 



a jV(<8*)/s*0 D 



(B.3) 



where S* denotes the isotopomer S in some given energy level 
from which the transition in question occurs, N(£*) refers to 
the number density of isotopomers £ in the excited energy level 
and Pb* refers to the oscillator strength or probability of the 
particular transition from the excited energy level. Now, at any 
given wavelength width A/l in Eq. IB.1I we see that 
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Now, ass uming LTE, we know from statistical mechanics 
dMihalaslll978h that a Boltzmann distribution of energy levels 
implies that where many instances of some microscopic system 
T occur, the population of a particular excited energy level T* 
is given by 



N(T*) = 



NCQgr* exp(^) 
Z(T) 



(B.5) 



where gy* is the statistical weight of the excited energy level, 
Ey* is its energy and Z(T) is the partition function of the sys- 
tem. Hence, 



Krg* a 



N(S) gs * exp(^)/ 2 



(B.6) 



Now, the number of molecules is related to the number density 
of the molecule's constituent atoms through the well-known 
Guldberg-Waage law, so 



N(S) 



N( b X)N( c Y) 
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where x(S) is the equilibrium constant of the f ormation- 
dissociation reacti on of the isotopomer £>, given JSchaded 
ll964lTaturrJf9 66) by 



Z(»X)Z( f Y)/4exp(ffi) 
Z(S) 



(B.8) 



where and are the dissociation energy and reduced mass 
respectively of isotopomer S. Hence, combining Eas. IB.6HB.8l 
we now see that 



N( b X)N( c Y) gfs* exp(^) 
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where gg+ and /g* have been combined into gf@*, the gf value 
of the line. 

So, the opacity scalefactor g K for an isotopomer S, where 
the reference isotopomer is J{ is then 



Sk = 



KB* 



Z(''X)Z('Y)^ 



N(rX)N('Y)gf^ exp( ggj ) 



iV(»X)g/ g »ZC'X)/ij, exp(^g^) _ 
iV(«X)g/«*Z(*X)4 exp(^al) 



y m :A T it 



(B.10) 



Now, we make the approximations g/g. w gfj\*, D<g ^ Dj\ and 
£3* » since and & are very nearly the same molecule 
so their transitions and energy levels will be almost identical. 
Further, assuming Z( b X) m Z(°X) due to the first-order depen- 
dence of atomic orbitals upon charge (and not mass), we arrive 
at 



(B.ll) 
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where the final scalefactor is independent of transition. The re- 
duced mass opacity corrections used for 13 C 16 and 12 C 18 
were therefore 0.9346 and 0.9293 respectively, wh e re nu clear 
masses were again sourced from lAudi A Wansti 3 1199.4 . In 
practice, the ratio of the square root terms enters the opacity 
anyway through the mass scalefactor, and the microturbulent 
velocity is not used in 3D (i.e. £ t = was used) so these can be 
disregarded also. The total opacity scalefactor applied in each 
case was therefore the product of the appropriate one of the re- 
duced mass correction factors and the ratio of reference ( 12 C or 
16 0) isotope to that involved in the line being modelled ( 13 C or 
18 0), i.e. 

N( b X) (K*^ _/*XUtf* 3 

Had the opacity scalefactor reflected abundance differences 
only, not intrinsic atomic (i.e. reduced mass) corrections, our 
adopted 3D-based ratios would have been I2 C/ 13 C = 92.9+^ 
and 16 0/ 18 = 516! 32 . 

The opacity scalefactors were then iteratively altered in cal- 
culations of the 13 C 16 and 12 C ls O line lists to reflect different 
isotopic ratios, in the same manner as the bulk carbon abun- 
dance was iterated to determine the ' 12 C abundance using the 
three sets of 12 C I6 lines. The same bulk carbon abundance as 
used in the final iteration of the weak 12 C 16 C calculation was 
used in each case, and the isotopic ratios after some iteration i 
given by 

^"X\ / n X\ io A ' 2c '™» 1 

' (B.13) 



b Xj; I'Xjj., 10^ 



Here Aiso, and A 12 Cfi na i are the abundance corrections pro- 
duced by the jth isotopomer iteration and the final iteration of 



whichever 12 C 16 list is used to indicate the 12 C abundance, 
respectively. 

B.3. Fractional Scalefactor 

Following opacity scalefactor convergence, a further iteration 
was also performed on each of the sets of 12 C 16 lines using 
a fractional scalefactor. This accounts for the fact that not all 
carbon in CO lines exists in 12 C 16 0. This scalefactor was cal- 
culated from the just derived 12 C/ 13 C and I6 0/ 18 ratios, such 
that 

N( 12 C 16 0) 



N(CO) ' 



and since 

N(CO) = N( 12 C 16 0) + N( 13 C 16 0) + N( 12 C I8 0) + N( 12 C 17 0) 
we see that 



Af( 1J C lb O) Af( lz C 18 Q) 
sW_l + N( u C 16 0) + JV( 12 C 16 0) + 

Af( 12 C 17 Q) 
+ A?( 12 C 16 0) 
= (1 + 13 C/ 12 C + 18 0/ 16 + 17 0/ 16 0) _1 . (B.14) 

Therefore, gf mc represents the fraction of the bulk carbon abun- 
dance actually indicated by 12 C I6 lines. Seeing as the 12 C 17 
lines are so weak in the ATMOS spectrum, accurate derivation 
of the 16 0/ 17 ratio would not have been possible in this study, 
so the value of 16 0/ 17 used in the calculation of fractional 
scale factors was the terrestrial value of ~2630 JCoplen et alJ 
l2002h . Being so incredibly small, the contribution of any dif- 
ference in the n O abundance between Earth and the Sun would 
have had almost no effect upon the resultant scalefactor. Using 
these fractional scalefactors, the final iterations of the 12 C 16 
line lists indicated bulk solar carbon abundances. 



Appendix C: Line Lists 

Refer to Tables lcTtiOl 
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Table C.l. Strong 12 C lfi O line list used for the Sect. |5] 
study, consisting of 31 lines, log gf data calculated from 
iGoorvitch & Chackerianldl994 . 



Table C.2. Weak 12 C I6 line list used for the Sect. © 
study, consisting of 13 lines, log gf data calculated from 
IGoorvitch & Chackerianl ( ll994 . 



Trans. Wavelength 
vac. (nm) 



IT 

(cm-') 



Exc. 
(eV) 



Trans. 



Wavelength 
vac. (nm) 



a 

(cm- 1 ) 



Exc. 
(eV) 



W A log a2 C 
(pm) (3D) 



1-0 P 6 
1-0 P 7 
1-0 R31 
1-0 R42 
1-0 R58 
1-0 R59 
1-0 R60 

1- R61 

2- 1 P 2 
2-1 R0 
2-1 R18 
2-1 R30 
2-1 R36 

2- 1 R41 

3- 2 R14 

3- 2 R27 

4- 3 R15 
4-3 R32 

4- 3 R38 
5^1 R19 
5^1 R28 

5- 4 R31 
5-4 R34 
5-4 R39 

5- 4 R40 
5^1 R48 

6- 5 R31 
6-5 R58 
6-5 R61 

6- 5 R68 

7- 6 R49 



4717.6910 
4726.7267 
4450.3204 
4398.4026 
4341.8489 
4339.0480 
4336.3331 
4333.7041 
4741.2775 
4715.7219 
4582.1879 
4510.9808 
4480.4509 
4457.5278 
4666.9654 
4584.3568 
4719.1825 
4614.8044 
4584.7513 
4752.1533 
4696.2029 
4679.3794 
4663.4526 
4638.8782 
4634.2564 
4600.7586 
4739.9862 
4627.2312 
4619.1075 
4603.5454 
4718.1039 



2119.681 
2115.629 
2247.029 
2273.553 
2303.166 
2304.653 
2306.096 
2307.495 
2109.136 
2120.566 
2182.364 
2216.813 
2231.918 
2243.396 
2142.720 
2181.331 
2119.011 
2166.939 
2181.143 
2104.309 
2129.380 
2137.036 
2144.334 
2155.694 
2157.844 
2173.555 
2109.711 
2161.120 
2164.920 
2172.239 
2119.496 



-4.171 
-4.104 
-3.412 
-3.277 
-3.129 
-3.121 
-3.113 
-3.106 
-4.350 
-4.648 
-3.355 
-3.134 
-3.053 
-2.996 
-3.292 
-3.011 
-3.146 
-2.818 
-2.742 
-2.955 
-2.788 
-2.745 
-2.703 
-2.642 
-2.631 
-2.550 
-2.674 
-2.394 
-2.372 
-2.322 
-2.411 



0.010 
0.013 
0.236 
0.428 
0.807 
0.834 
0.862 
0.891 
0.267 
0.266 
0.346 
0.485 
0.579 
0.670 
0.577 
0.705 
0.843 
1.031 
1.129 
1.130 
1.229 
1.270 
1.316 
1.400 
1.418 
1.579 
1.521 
2.066 
2.146 
2.347 
2.093 
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1- 
P 96 
P 98 
R 98 
R 99 
R100 
R101 
R103 
R105 
R107 
R109 

2- 1 
P 97 
P 98 
P 101 



6112. 
6162. 
4297. 
4297. 
4298. 
4299. 
4301. 
4303. 
4306. 
4309. 



9978 
3540 
1649 
8564 
6396 
5151 
5437 
9460 
7253 
8852 



6225.7984 
6251.2034 
6329.1932 



1635.859 
1622.757 
2327.116 
2326.741 
2326.317 
2325.844 
2324.747 
2323.449 
2321.950 
2320.247 

1606.220 
1599.692 
1579.980 



-3.053 
-3.046 
-2.889 
-2.883 
-2.879 
-2.873 
-2.867 
-2.857 
-2.848 
-2.842 



2.154 
2.242 
2.242 
2.286 
2.331 
2.376 
2.467 
2.560 
2.654 
2.750 



7.74 
6.72 
8.46 
7.72 
7.21 
6.69 
5.74 
4.95 
4.02 
3.38 



-2.755 2.443 7.67 
-2.752 2.486 7.11 
-2.742 2.619 5.61 



8.37 
8.38 
8.41 
8.39 
8.40 
8.39 
8.40 
8.40 
8.39 
8.40 

8.39 
8.39 
8.38 



Table C.3. 



13 C le O line list used for the Sect. |6] study, 
of 16 lines, log gf data calculated from 



consisting 

IGoorvitch & Chackeriatj ( ll994 . 



Trans. 



Wavelength 
vac. (nm) 



a 

(cm" 1 ) 



Exc. 
(eV) 



W.i log ei3 C 
(pm) (3D) 



1- 
R 42 
R 43 
R 44 
R 47 
R 55 
R 56 
R 54 

2- 1 
R 46 
R 47 
R 49 
R 53 
R 62 
R 67 
R 77 
R 79 

3- 2 
R 62 



4501.6474 
4497.3756 
4493.1896 
4481.1460 
4452.7709 
4449.6042 
4456.0218 



4539. 
4535. 
4528. 
4513. 
4486. 
4474. 
4457. 
4454. 



7468 
7888 
1342 
8671 
8117 
7857 
1692 
6795 



2221.409 
2223.519 
2225.591 
2231.572 
2245.793 
2247.391 
2244.154 

2202.766 
2204.688 
2208.415 
2215.395 
2228.754 
2234.744 
2243.577 
2244.830 



-3.297 
-3.286 
-3.275 
-3.245 
-3.173 
-3.165 
-3.182 



0.409 
0.429 
0.448 
0.511 
0.695 
0.720 
0.671 



-2.963 0.745 
-2.951 0.766 
-2.932 0.809 
-2.896 0.901 
-2.827 1.132 
-2.790 1.275 
-2.726 1.591 
-2.714 1.660 



4.26 
4.25 
4.24 
4.02 
2.91 
2.77 
2.78 

3.92 
4.11 
3.79 
3.46 
2.23 
1.80 
0.99 
0.89 



6.43 
6.44 
6.45 
6.46 
6.43 
6.42 
6.39 

6.42 
6.46 
6.45 
6.47 
6.43 
6.45 
6.46 
6.48 



4542.4821 2201.440 -2.658 1.380 1.88 6.46 
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Table C.4. 12 C ls O line list used for the Sect. |5] study, 
consisting of 15 line s. log gf data calculated from 



iGoorvitch & Chackerianld 199-4 



Trans. 


Wavelength 


cr 




vac. (nm) 


(cm- 1 ) 


1-0 






P 26 


5033.4006 


1986.728 


P 27 


5044.7281 


1982.267 


P 34 


5127.5809 


1950.237 


P 43 


5243.7159 


1907.045 


R 48 


4486.1030 


2229.106 


R 49 


4482.3365 


2230.979 


R 50 


4478.6546 


2232.813 


2-1 






P 25 


5084.5153 


1966.756 


P 29 


5130.8774 


1948.984 


P 31 


5154.8355 


1939.926 


3-2 






P 35 


5270.2261 


1897.452 


P 40 


5336.4145 


1873.917 


4-3 






P 21 


5167.4916 


1935.175 


P 23 


5190.1890 


1926.712 


5-4 






P 34 


5393.2644 


1854.165 



log gf 



Exc. 
(eV) 



log eis 
(pm) (3D) 



3.572 


0.159 


1.09 


6.01 


3.556 


0.171 


1.02 


5.97 


3.462 


0.269 


0.89 


5.92 


3.370 


0.427 


0.88 


6.00 


3.237 


0.530 


0.74 


5.94 


3.228 


0.552 


0.74 


5.95 


3.219 


0.574 


0.72 


5.95 


3.293 


0.405 


1.06 


5.99 


3.233 


0.455 


1.00 


5.96 


3.205 


0.482 


0.98 


5.95 


2.987 


0.796 


0.78 


5.97 


2.932 


0.879 


0.80 


6.02 


3.077 


0.871 


0.45 


5.91 


3.039 


0.891 


0.59 


6.01 


2.788 


1.278 


0.38 


5.99 



Table C.5. Low excitation (LE) 12 C le O line list used for the 
Sect.EQstudy, consisting of 15 lin es, log gf data calculated from 
IGoorvitch & Chackerianl (11 994h. 

Trans. Wavelength cr logg/ Exc. W,i logei2 C 
vac. (nm) (cm 4 ) (eV) (pm) (3D) 



2- 


-1 














P 


2 


4741.2775 


2109.136 


-4.350 


0.267 


20.79 


8.55 


R 





4715.7219 


2120.566 


-4.648 


0.266 


15.14 


8.47 


3- 

P 

4 


-2 
2 
-3 


4801.2476 


2082.792 


-4.181 


0.530 


18.37 


8.51 


P 


2 


4862.5943 


2056.515 


-4.063 


0.789 


15.37 


8.47 


P 


3 


4871.6187 


2052.706 


-3.886 


0.790 


18.37 


8.49 


P 


4 


4880.7591 


2048.862 


-3.762 


0.792 


20.84 


8.54 


R 





4836.2087 


2067.735 


-4.362 


0.787 


10.53 


8.43 


R 


1 


4827.6411 


2071.405 


-4.060 


0.788 


15.50 


8.48 


5-4 














P 


1 


4916.3051 


2034.048 


-4.273 


1.044 


7.86 


8.42 


P 


3 


4934.5350 


2026.533 


-3.799 


1.046 


14.68 


8.44 


P 


4 


4943.8275 


2022.724 


-3.674 


1.048 


17.11 


8.48 


P 


6 


4962.7720 


2015.003 


-3.499 


1.053 


20.05 


8.50 


R 





4898.5440 


2041.423 


-4.272 


1.043 


7.58 


8.39 


R 


2 


4881.2454 


2048.658 


-3.793 


1.045 


14.87 


8.45 


R 


4 


4864.4047 


2055.750 


-3.569 


1.048 


18.93 


8.50 
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Tabl e C.6. First overtone (Av = 2) 12 C 16 line list used for the Sect.[5]study, consisting of 66 lines, log gf data calculated from 
iGoorvitch & Chackerianldl994 . 



Trans. Wavelength a logg/ Exc. W A logei2 C 
vac. (nm) (cm" 1 ) (eV) (pm) (3D) 



2 


-0 














P 


7 


2363.1249 


4231.685 


-6.223 


0.013 


1.04 


8.34 


P 


9 


2368.0105 


4222.954 


-6.119 


0.021 


1.31 


8.35 


P 


13 


2378.3082 


4204.670 


-5.971 


0.043 


1.73 


8.36 


P 


20 


2398.0498 


4170.055 


-5.799 


0.100 


2.13 


8.35 


P 


21 


2401.0522 


4164.841 


-5.780 


0.110 


2.22 


8.36 


P 


22 


2404.1007 


4159.560 


-5.762 


0.120 


2.23 


8.35 


P 


32 


2437.1831 


4103.098 


-5.622 


0.251 


2.22 


8.35 


P 


38 


2459.3740 


4066.075 


-5.559 


0.352 


2.00 


8.35 


P 


41 


2471.1542 


4046.692 


-5.532 


0.408 


1.85 


8.34 


P 


42 


2475.1846 


4040.103 


-5.523 


0.428 


1.84 


8.35 


P 


46 


2491.8344 


4013.108 


-5.491 


0.512 


1.61 


8.35 


P 


53 


2523.0648 


3963.434 


-5.443 


0.676 


1.24 


8.36 


P 


54 


2527.7505 


3956.087 


-5.435 


0.701 


1.13 


8.33 


R 


9 


2328.4851 


4294.638 


-6.024 


0.021 


1.63 


8.36 


R 


13 


2322.0848 


4306.475 


-5.867 


0.043 


2.11 


8.36 


R 


16 


2317.7138 


4314.597 


-5.775 


0.065 


2.47 


8.37 


R 


19 


2313.7084 


4322.066 


-5.695 


0.090 


2.68 


8.36 


R 22 


2310.0671 


4328.879 


-5.625 


0.120 


2.91 


8.36 


R 23 


2308.9340 


4331.003 


-5.606 


0.131 


2.94 


8.36 


R 25 


2306.7886 


4335.031 


-5.564 


0.155 


3.00 


8.36 


R 26 


2305.7763 


4336.934 


-5.545 


0.167 


3.08 


8.37 


R 41 


2295.4058 


4356.528 


-5.301 


0.408 


2.86 


8.36 


R 65 


2297.7720 


4352.042 


-5.041 


1.009 


1.22 


8.34 


R 66 


2298.3892 


4350.873 


-5.032 


1.039 


1.25 


8.37 


P 


— J 

13 


2408.3440 


4152.231 


-5.487 


0.309 


2.61 


8.36 


P 


15 


2413.8559 


4142.749 


-5.430 


0.322 


2.83 


8.36 


P 


16 


2416.6808 


4137.907 


-5.403 


0.330 


2.93 


8.36 


P 


18 


2422.4694 


4128.019 


-5.357 


0.346 


3.13 


8.37 


P 


22 


2434.6085 


4107.437 


-5.279 


0.385 


3.30 


8.36 


P 


23 


2437.7618 


4102.124 


-5.262 


0.396 


3.30 


8.35 


P 


29 


2457.6960 


4068.851 


-5.174 


0.471 


3.38 


8.36 


P 


30 


2461.1901 


4063.075 


-5.162 


0.485 


3.42 


8.36 



Trans. 


Wavelength 


cr 


log,?/ 


Exc. 


W A 


log ei2 C 






vac. (nm) 


(cnr 1 ) 




(eV) 


(pm) 


(3D) 


P 


32 


2468.3280 


4051.326 


-5.138 


0.514 


3.31 


8.35 


P 


33 


2471.9723 


4045.353 


-5.126 


0.530 


3.30 


8.36 


P 


36 


2483.2111 


4027.044 


-5.095 


0.579 


3.11 


8.35 


P 


38 


2490.9617 


4014.514 


-5.075 


0.614 


3.01 


8.35 


P 


40 


2498.9218 


4001.726 


-5.057 


0.651 


2.90 


8.35 


P 


58 


2580.5013 


3875.216 


-4.928 


1.065 


1.51 


8.35 


P 


60 


2590.7333 


3859.911 


-4.917 


1.120 


1.36 


8.35 


R 


8 


2359.4419 


4238.290 


-5.590 


0.283 


2.14 


8.35 


R 


13 


2351.2311 


4253.091 


-5.385 


0.309 


3.10 


8.36 


R 


14 


2349.7148 


4255.836 


-5.353 


0.315 


3.26 


8.36 


R 


15 


2348.2403 


4258.508 


-5.322 


0.322 


3.37 


8.35 


R 


16 


2346.8074 


4261.108 


-5.293 


0.330 


3.50 


8.35 


R 


17 


2345.4162 


4263.636 


-5.266 


0.338 


3.69 


8.36 


R 


18 


2344.0666 


4266.090 


-5.234 


0.346 


3.79 


8.36 


R 


19 


2342.7584 


4268.473 


-5.214 


0.355 


3.90 


8.36 


R 


20 


2341.4918 


4270.782 


-5.190 


0.365 


4.03 


8.36 


R 


22 


2339.0825 


4275.181 


-5.146 


0.385 


4.21 


8.37 


R 


23 


2337.9399 


4277.270 


-5.124 


0.396 


4.25 


8.36 


R 


24 


2336.8385 


4279.286 


-5.104 


0.407 


4.37 


8.37 


R 


31 


2330.2818 


4291.327 


-4.975 


0.499 


4.89 


8.41 


R 


32 


2329.5096 


4292.749 


-4.963 


0.514 


4.61 


8.38 


R 


35 


2327.4396 


4296.567 


-4.914 


0.562 


4.50 


8.37 


R 


37 


2326.2652 


4298.736 


-4.886 


0.596 


4.46 


8.37 


R 


38 


2325.7397 


4299.707 


-4.870 


0.614 


4.54 


8.39 


R 


40 


2324.8122 


4301.423 


-4.845 


0.651 


4.32 


8.37 


R 


41 


2324.4103 


4302.167 


-4.830 


0.670 


4.28 


8.37 


R 


56 


2323.3645 


4304.103 


-4.654 


1.012 


2.84 


8.36 


R 


63 


2326.1201 


4299.004 


-4.583 


1.206 


2.19 


8.37 


R 


64 


2326.6855 


4297.959 


-4.574 


1.235 


2.08 


8.37 


R 


65 


2327.2943 


4296.835 


-4.564 


1.265 


2.14 


8.40 


R 


68 


2329.3824 


4292.984 


-4.535 


1.357 


1.83 


8.40 


R 


70 


2330.9940 


4290.015 


-4.517 


1.421 


1.58 


8.38 


R 


77 


2338.0411 


4277.085 


-4.456 


1.657 


0.99 


8.36 


R 


80 


2341.7430 


4270.324 


-4.431 


1.765 


0.79 


8.35 
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